Pnictides as frustrated quantum antiferromagnet close to a quantum phase transition 



Gotz S. Uhrig *It| Michael Holt, Jaan Oitmaa, and Oleg P. Sushkov 
School of Physics, University of New South Wales, Kensington 2052, Sydney NSW, Australia 



o 
o 

(N 
D 



Rajiv R.P. Singh 

University of California Davis, CA 95616, USA 
(Dated: February 9, 2009) 

We present a detailed description of the dynamics of the magnetic modes in the recently discov- 
ered superconducting pnictides using reliable self-consistent spin- wave theory and series expansion. 
Contrary to linear spin-wave theory, no gapless mode occurs at the Neel wave vector. We discuss the 
scenario that the static magnetic moment is strongly reduced by magnetic fluctuations arising from 
the vicinity to a quantum phase transition. Smoking gun experiments to verify this scenario are 
proposed and possible results are predicted. Intriguingly in this scenario, the structural transition 
at finite temperature would be driven by an Ising transition in directional degrees of freedom. 
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The seminal discovery of superconductivity in the pnic- 
tide family of materials at relatively high temperatures[l] 
has led to tremendous excitement and research activity. 
This discovery raises many fundamental questions. Fore- 
most among them is whether high temperature supercon- 
ductivity in these materials is in some fundamental sense 
closely related to those in the cuprate family of materi- 
als. Indeed the quasi-two dimensional layered structure 
for the two families and antiferromagnetism in the par- 
ent compounds suggests potential similarities. However, 
many doubts have also been raised about any correspon- 
dence, such as: Is the origin and nature of spin fluctu- 
ations in the two families related given that the parent 
compounds are metallic in the pnictides whereas they are 
insulating in the cuprates? Are the pnictide materials 
even strongly correlated or are Local Density Approx- 
imation based approaches adequate? Are spin models 
appropriate for describing spin-fluctuations in these ma- 
terials? In order to address these very basic questions, 
it is necessary to have detailed quantitative comparisons 
between theory and experiments. 

On the experimental front, magnetic long range order 
was established in LaFeAsOi_a;Fa; by neutron scattering 
(NS) [2| and by muon spin resonance (/zSR) The NS 
provides evidence for a columnar antiferromagnetic or- 
dering with a staggered magnetic moment of 0.36(5)/iB- 
For simplicity we consider here only the square lattice 
which is formed by the Fe ions ignoring a small or- 
thorhombic and even monoclinic structural distortion. 
Along the a axis the spin directions alternate whereas 
they are the same along the b axis, see Fig. [Ija). The 
/iSR also provides evidence that the spin order is com- 
mensurate but with a small staggered moment of 0.25/iB- 
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First results on the dispersion of the magnetic excitations 
have just become available [1, A tiny anisotropy gap 
is found to be « 6meV, the spin wave velocity v±_ per- 
pendicular to the stripes to be 205±20meV in units of 
l/g, where g is the inverse Fe-Fe distance, and a small 
interplane coupling Jz is found to be « 5meV. Results 
for the parallel spin wave velocity are not available so 
far, but they are expected soon. 
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FIG. 1: (color online) Panel (a) considered spin pattern; panel 
(b) staggered magnetization for S = 1/2 and 5 = 1 as func- 
tion of the ratio of the couplings. For comparison the re- 
sults of linear spin wave theory (LSW) are included. The 
self-consistent spin wave theory is nicely supported by series 
expansion about the Ising limit (symbols). 

Theoretically, the columnar antiferromagnetic ordering 
was also found to be the most stable in band structure 
calculation @, So there is agreement on the static 
structure. But the smallness of the staggered magnetic 
moment is a matter of controversy. On the one hand, 
band structure results indicate a local moment of up 
to 2.3^B 0, H, Hi- This has led to the suggestion that 
the magnetic fluctuations themselves strongly reduce the 
static local moment [13] • We will show that this scenario 
implies that the pnictides are in the direct vicinity of a 
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quantum phase transition. On the other hand, there are 
studies suggesting that the strong reduction of the local 
magnetic moment can be explained by electronic effects 
like hybridization, spin-orbit coupling, and a particular 
low symmetry [111 Il2| . Hence, there are two different 
scenarios: (i) the local static moment is reduced by the 
magnetic fluctuations. In this case the ratio of couplings 
must be fixed to an appropriate value, (ii) the local elec- 
tronic orbitals account for the sizable reduction so that 
the magnetic couplings are not determined by the value 
of the magnetic moment. 

Together with upcoming experimental results, our 
work will help decide on the degree of strong correlation 
and on the closeness to a quantum phase transition. 

The quantitative goal of the present work is three- 
fold. First, we provide a theory for the magnetic ex- 
citations based on a minimal spin model, namely the 
J1-J2 Heisenberg model with spin S" = 1 (5 = 1/2 re- 
sults are shown for comparison and to justify our approx- 
imations). Measurements of the spin wave dispersion up 
to its maxima will further support the spin-Hamiltonian 
approach whose validity in turn shows that strong cor- 
relations dominate the pnictides. We show that a siz- 
able finite energy of the spin waves is to be expected at 
q = (1,1) where we denote all wave vectors in units of 
7r/(/. This is in stark contrast to the results of linear spin- 
wave theory (LSW) where a vanishing spin wave energy is 
predicted [l3,[l3]- Second, we discuss the possibility that 
the magnetic fluctuations reduce the static local moment. 
Third, we make quantitative predictions of the dispersion 
and of the anisotropy in spin wave velocities along (w||) 
and across (wj_) the magnetic stripes. Measurements of 
v\\/v± can be used to determine the magnetic frustration 
J2/J1. The spin- wave spectra over the full Brillouin zone 
(BZ) show clear differences between a system deep in the 
columnar phase and one close to a quantum phase tran- 
sition, where magnetic fluctuations dramatically reduce 
the static moment. This provides a robust experimental 
way to distinguish the two scenarios (i) and (ii). 

The parent compound LaFeAsO is not a Mott- 
Hubbard insulator. It is rather a bad metal or semimetal 
without a Drude peak in the conductivity 0, [l3| • Even 
in the magnetically ordered phase the entire Fermi sur- 
face is not gapped @, Q. Still, in the undoped system 
the magnetic excitations are long-lived as they appear as 
sharp peaks [3] even at high energies [s^. So it is justi- 
fied, though not undisputed , to start with a model of 
well-localized spins coupled by Heisenberg interactions 
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where («, j) stands for nearest neighbors (NN) and 
for next-nearest neighbors (NNN). The spin operators 
read and represent 5 = 1. This is not completely 
obvious in view of the complicated local electronic sit- 
uations. But the success of two-band models strongly 



indicates that S* = 1 is the appropriate choic e [itI . [l^ . 
Furthermore, the band structure results 0, d, Qwith a 
local moment of up to 2.3/xb also indicates that there can 
be up to two electrons aligned. While our choice ([T]) ne- 
glects life time effects due to the decay into particle-hole 
pairs (Landau damping) we expect that the collective 
magnetic excitations and their dispersion are captured. 

The choice of a J1-J2 Heisenberg model is justified be- 
cause the superexchange is realized mostly via the As ion 
which sits in the middle of each Fe plaquette. Hence a NN 
contribution and a NNN contribution is to be expected 
[lot . Indeed, band structure calculations show that the 
NN and the NNN coupling are sizable 0, H [ll- Both 
Ji and J2 turn out to be antiferromagnetic, i.e., positive 
and very similar in value. For this reason, we will choose 
for scenario (ii) the ratio J1/J2 = 1. 

The technique employed is self-consistent spin wave 
theory. It has been shown previously that this approxi- 
mation works extremely well in the columnar, stripe-like 
phase Him for 5=1/2, see also Figs.IHb) andH^b). 
Because spin wave theory can be derived as a 1/5 expan- 
sion, the results should only improve for 5=1. 

We have used the Dyson-Maleev as well as the 
Schwinger boson representation 22| . Both yield the same 
result on the level of self-consistent mean-field theory. In 
the symmetry broken phase the dispersion reads 



w(q) = CJ2\/A^-B^ (2a) 
A = A + x/3cos(7rgb) (2b) 
B — 2cos{TTqa) cos{Trqb) + xacos{TTqa) (2c) 

with X — J1/J2. The expectation values C, a and /3 are 
determined from the self-consistency conditions 

25 = dEMF/d{CJ2X) + 2m (3a) 

2CJ2 = -dEMF/d{CJ2)+4J2m (3b) 

CaJi = -dEMF/d{Ca) + 2Jim (3c) 

CpJi = dEMF/d{CP) + 2Jim, (3d) 

where m is the staggered magnetization and the mean- 
field energy per spin 



E, 



MF 



(27r) 



(c^(q) - ACJ2)d^q 



(4) 



BZ 



is used. One integration can be done analytically, the 
second numerically. The self-consistency is solved by it- 
eration. Gaplessness at zero wave vector implies A = 
2 -I- x{a — f3) so that ([3a|) is used to determine m. 

The resulting magnetizations are shown in Fig. [ijb). 
Note the extremely fast vanishing of the magnetization 
if X approaches 2. The vicinity of the classical first order 
instability at a; = 2 [23| makes m{x) resemble a square 
root a,s X ^ Xc where m vanishes {xc\s=i/2 = 1-8057 and 
Xc\s=i — 1-9836) within the approximation. One may 
speculate that this is due to the Ising-type transition re- 
lated to the breaking of directional symmetry 1^ 24 , 2^ , 
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but so far we cannot draw a definitive conclusion on this 
point. But in the light of the structural phase transition 
occurring before 0,11] or at [H, [S^l the magnetic phase 
transition this aspect is experimentally very interesting. 
The structural transition could easily be driven by the 
Ising transition in the directional degrees of freedom. 

For completeness, we also include results for negative 
Ji. Around x ~ —2 another instability is expected [11] ■ 
But interestingly it does not lead to any precursors in 
the sublattice magnetization as found from self-consistent 
spin-wave theory and series expansion. We attribute this 
behavior to the fact that the classical instability to the 
ferromagnetic phase is completely first order in the sense 
that there are no precursive fluctuations because the fer- 
romagnetic phase is free from quantum fluctuations. 
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FIG. 2: (color online) Panel (a) ratio of the spin wave veloci- 
ties parallel and perpendicular to the spin stripes, i.e., parallel 
is in 6 direction and perpendicular is in a direction, see Fig. 
[TJa). We stress that , though small, does not vanish where 
the magnetization vanishes. For comparison, LSW data is 
also included. Panel (b) depicts the spin wave energy aj(q) 
at q = (1, 1) and q = (0, 1), respectively, in units of n/g 
{g lattice constant of the assumed square lattice). Note that 
i^{<l) , though small, does not vanish where the magnetization 
vanishes. The results are nicely corroborated by exemplary 
series expansion data (symbols). 

The breaking of the directional symmetry implies that 
the spin wave velocities depend on direction, see Fig. 
[2)Ja). This quantity is a much more robust probe for 
the value of the ratio x = J1/J2. The magnetization 
depends on matrix elements which in turn can depend 
on itinerancy, hybridization and other effects. Energies 
in contrast only depend on the Hamiltonian and thus are 
much less ambiguous. 

Furthermore, one notes that ^ implies that there is a 
finite excitation energy at q — (±1. ±1) which is equiv- 
alent to q = (0, ±1) if x{a — /3) > 0. The results are 
plotted in Fig. ^h) . We stress that no dependence of the 
bare coupling Ji on the bond direction is required. 

Now we turn to the two scenarios presented. In sce- 
nario (i) we attribute the reduction of m to renormaliza- 
tion by the spin fluctuations. Equating the ratio 0.36/2.3 
as from experiment |3| and band structure theory [3, 0| 
to m in Fig. Wljo) leads to Ji = 1.978J2. From Fig. [ijjb) 
it is obvious that considerable fine-tuning is needed. The 
resulting dispersion is shown in the upper panel in Fig. 



FIG. 3: Upper panel: dispersion of scenario (i) with S — 1, 
J2 — 34meV and J1/J2 = 1.978; lower panel: dispersion of 
scenario (ii) with S = 1, J2 = SSmeV and J1/J2 = 1. J2 is 
chosen such that v± equals the experimental value Wave 
vectors in units of vr/g. 



[21 The lower panel shows the dispersion for Ji = J2 
as suggested by band structure calculations J3j \^ ■ In 
both cases the overall scale of the coupling is adjusted to 
fit to the measured spin velocity v± (J]. We stress that 
the value for J2 of about 33meV agree very well with the 

Em. 



estimates from band structure calculations 

Comparing the two panels of Fig. [31 the difference in 
the dispersion of the spin waves in both scenarios is strik- 
ing. Thus a measurement of the energetically higher ly- 
ing modes will easily distinguish both scenarios. To fa- 
cilitate the distinction we plot in Fig. [H the dispersion 
along a generic path in the BZ. The perpendicular spin 
wave velocity v± is fixed to its experimental value Q 
by the appropriate choice of J2. Obviously, the motion 
across the stripes (along a*) is the same for both sets of 
parameters. The important difference occurs in the mo- 
tion along the stripes (along b*). If the spin fluctuations 
renormalize the magnetic moment the spin mode along 
b* is extremely soft. Interestingly, this softness opens 
an additional channel for the magnons to decay at ener- 
gies of about twice the energy at (0, 1), i.e., above about 
40-50 meV. So in this scenario significant line widths in 
inelastic neutron scattering are to be expected. 
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FIG. 4: (color online) Dispersions of both scenarios for a 
generic path through the BZ. They coincide along (0, 0) 
(1,0), but differ strongly along (0,1) (0,0). In LSW the 
dispersions at (0, 1) would vanish spuriously. Wave vectors 
are given in units oi n/g. 

In conclusion, we presented a quantitative theory for 
the dispersion of the spin waves in the recently discovered 
superconducting pnictides. It is based on self-consistent 
spin wave theory and series expansion for the S — 1 Ji- 
J2 Heisenberg model. Measurements at higher energies 
will further support a model of localized spins, for first 
evidence see Refs. |^ [2^ We predict a strong anisotropy 
of the spin wave velocities and a finite excitation energy 
at the wave vectors (0, 1) and (1, 1). 

Two scenarios for the strong reduction of the local 
magnetic moment are considered. The scenario (i) at- 
tributes the reduction to the magnetic fluctuations. We 
point out that the strongly varying static moments from 
Q.2^^lB [i over QM^b 01 to C.g^ie ^ and 0.9/xb 
finds a natural explanation if the pnictides are close to 
the quantum phase transition at a; w 2 where the renor- 
malized magnetization changes very rapidly on small pa- 
rameter changes. This scenario implies the fascinating 
aspect that the pnictides realize a spin-isotropic system 
which displays an Ising transition in the orientation of 
its ferromagnetic stripes [3, 111, [111 ■ 

The alternative scenario (ii) attributes the low mag- 
netic moment to the local electronic configuration. Then 
J\l J2 ~ 1 is plausible which does not renormalize the 
magnetic values sizably (about 18%). 

The measured anisotropics will allow one to decide how 
close the system is to a quantum phase transition with 
fascinating features such as directional Ising transitions. 
A quantitative understanding of spin fluctuations in the 
parent materials would help clarify one of the most in- 
triguing issues in the field, namely, the similarities and 
differences between the pnictide and the cuprate family 
of materials at a fundamental level. 
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